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The ground state of Sn/Si(lll) and Sn/Ge(lll) surface a-phases is reexamined theoretically, 
based on ab — initio calculations where correlations are approximately included through the orbital 
dependence of the Coulomb interaction (in the local density + Hubbard U approximation). The 
effect of correlations is to destabilize the vertical buckling in Sn/Ge(lll) and to make the surface 
magnetic, with a metal-insulator transition for both systems. This signals the onset of a stable 
narrow gap Mott-Hubbard insulating state, in agreement with very recent experiments. Antiferro- 
magnetic exchange is proposed to be responsible for the observed F-point photoemission intensity, 
as well as for the partial metallization observed above above 60 K in Sn/Si(lll). Extrinsic metal- 
lization of Sn/Si(lll) by, e.g. alkali doping, could lead to a novel 2D triangular superconducting 
state of this and similar surfaces. 

PACS numbers: 71.30.-l-h, 73.20.-r, 73.20.At, 75.70.Rf, 75.70.Ak, 71.27.+a,74.78.Na 



Metal-insulator transitions and strongly correlated 
states of electrons in the genuinely two dimensional (2D) 
"dangling bond" surface states of semiconductors have 
long been sought after [l| but seldom realized and char- 
acterized. On these surfaces, band physics usually takes 
over, driving large structural reconstructions, wich re- 
move metallicity with a strong energy gain [3]. The re- 
sulting passivation of surface states unfortunately also 
removes alternative and more interesting phases, includ- 
ing charge density waves (CDW) and spin density waves, 
Mott-Hubbard insulators (Mis), and possibly 2D super- 
conductivity 0]. 

Plummer and others [1, d, @, 0, S] however called 
attention to the a-phase surfaces, obtained by cover- 
ing an ideal metallic (111) semiconductor surface with a 
V3x V3i?30° (V3) triangular array of group IV adatoms. 
These systems, where adatoms lie as far apart as ~ 7 
A, possess a very narrow half-filled adatom surface state 
band, making them ideally prone to various instabilities 
and to strong correlations. A low temperature CDW-like 
reversible 3x3 periodic surface adatom distortion was 
indeed reported in metallic Pb/Ge(lll) and Sn/Ge(lll) 
0, 8] whereas undistorted triangular MI states appeared 
to prevail in isoelectronic surfaces like Si/SiC(0001) [7] 
andK:Si(lll)\/3-J5 Q. 

Density functional calculations in the local density ap- 
proximation (LDA) supported this diversity of behavior, 
indicating that a large ~ 0.3 A periodic "up-down" dis- 
tortion [lol . [lH akin to a valence disproportionation 
[l^ should be the ground state of Sn/Ge(lll), against 
an undistorted magnetic insulator prevailing in a large 
gap system like Si/SiC(0001) [HIm, 15]. Interestingly, 
the intermediate case of Sn/Si(lll) failed to fall neatly 
on either side of this divide. Within LDA, this surface 
is equally close to undistorted magnetism (evolving to 
an insulator as we shall see later) as it is to a 3 x 3 
distorted metal [l3|; but neither state is actually stable 



in the strain free surface As we will show presently, 
T = LDA seems paradoxically to describe better the be- 
havior of these surfaces at non-zero temperature, whereas 
a better account of correlations is needed to describe their 
actual ground state. 

Systems realizing in 2D a spontaneous transition be- 
tween these two types of state, namely the distorted (or 
undistorted) nonmagnetic band metal and the magnetic 
undistorted MI are very interesting to pursue. Among 
other things the latter constitute the building block 
of many important strongly correlated systems, includ- 
ing (in the square lattice version) cuprate superconduc- 
tors. Model studies do indicate that distorted metal- 
undistorted MI transitions are to be expected in the 
present system as a function of parameters 15|. Cor- 



relations in Sn/Si(lll) and Sn/Ge(lll) surfaces were 
discussed by Flores et al. [l6l |. who however concluded 
against a transition to a MI ground state. In fact, un- 
til recently no such transition was actually observed in 
ck-phases. 

Within the last few months that situation changed 
drastically. Cooling Sn/Ge(lll) below 20 K apparently 
turns it from a 3 x 3 distorted metal to an undistorted 
\/3 insulator, presumably a MI [I7i|. Equally striking, 
Sn/Si(lll) is now shown by Modesti's group [3] to 
turn continuously from an undistorted metal to a nar- 
row gap insulator, again presumably a MI, below 60 
K. If Sn/Ge(lll) is intriguing enough due to the disap- 
pearance of structural distortion at low T, the result on 
Sn/Si(lll) is no less puzzling. Low temperature photoe- 
mission [3) US] finds Sn-related surface bands below Ep 
at the F-point, as though folded over from the if -point. 
This is an intriguing but clear indication of 3 x 3 period- 
icity, for in the -s/S surface there are no such filled surface 
states at T [20| . Structural tools including scanning tun- 
neling microscopy (STM) and photeoelectron diffraction 
show only \/3 periodicity, so that the 3x3 motif is not 
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FIG. 1: Sn/Si(lll) (dashed) and Sn/Ge(lll) (solid) cal- 
culations for increasing correlation parameter U: (a) mini- 
mum band gap, (b) spin moment per adatom, (c) Sn verti- 
cal height measured over the outer semiconductor plane, (d) 
magnetic exchange splitting. Lines are guides for the eye. 



structural. For both surfaces therefore, the LDA predic- 
tions of metallic ground states turn out to be in error. 

Unraveling this situation calls for a renewed theoretical 
effort. Calculations should at the same time be of first 
principles quality, so as to permit total energy compar- 
isons, but also treat correlations more accurately than 
LDA, so as to identify Mis if and when present. The 
LDA+U approach, while still a mean field approxima- 
tion (thus for example replacing a MI with a fictitious 
magnetic band insulator) does satisfy these criteria 1211 
and is suitable to describe quantitatively surface Mis [22] . 
We therefore conducted a series of accurate LDA and 
LDA-fU calculations for Sn/Si(lll), and for Sn/Ge(lll). 

The geometries considered were periodic slabs consist- 
ing of three (111) semiconductor bilayers, H-saturated at 
the bottom, and with the 1/3 monolayer of Sn adsorbed 
at T4 sites on the top surface. We used the plane wave 
implementation of density functional theory [l^j HI] in 
the gradient corrected local spin density approximation 
(LSD A), extended to include a Hubbard U (LSDA+U). 
Beginning with an undistorted geometry, (one Sn per 
cell) we turned on an increasing electron-electron onsite 
repulsion U ranging from zero to an estimated full value 
of J7 ~ 4 eV[2J] for two electrons occupying the same Sn 
5pz orbital - this sing le orbital constituting about 50% of 
the surface state 11 1|. The effect of a finite U parameter 
in LSDA-I-U is to favor integer occupancy of this orbital, 
which is the crucial effect of strong correlations. 

FiglD summarizes the main physical quantities 
calculated at zero-temperature for Sn/Si(lll) and 
Sn/Ge(lll), illustratively shown for increasing U val- 
ues. Sn/Si(lll), metallic and weakly magnetic already 
at J7 ~ 0, develops a stronger magnetization and even- 
tually a metal insulator transition J7 ~ 2 eV, where the 
spin moment per adatom site saturates to l/is. The in- 



I 30 



:> 20 

-10 



1 


1 




' 1 


1 ' 1 












- 

- 






c 


• / 






,0 














1 


1 




1 


1 , 1 



-0.6 



0.4 



-0.2 





A (A) 



0.2 



0.4 



0.6 



FIG. 2: Sn/Ge(lll): Total energy per adatom as a function 
of 3 X 3 up-down distortion A. LDA (solid line) and LSDA-I-U, 
(7=4 eV (dashed line). The zero of energy is at A = 0.0 A 
in both cases. Inset: Schematic of \/3 x \/3R30° unit cell 
(dashed line) and of 3 x 3 unit cell (solid line). Black spheres 
represent Sn adatoms, grey spheres Si (Ge) atoms of the first 
and second layer. 



sulating gap reaches 0.3 eV at the final realistic value 
f7 ~ 4 eV, where the exchange splitting between oppo- 
site spins orientations becomes as large as 0.8 eV. The 
insulating state gains Coulomb energy, but loses band en- 
ergy relative to the metallic state. This is also reflected 
indirectly by the predicted downward geometrical relax- 
ation of the Sn adatoms by about 0.06 A closer to the Si 
substrate, upon going from metal at to insulator. Ver- 
tical adatom positions control hybridization between the 
5pz dangling bond orbital and the underlying Si-Si an- 
tibonding state [HI, [2^ . That hybridization is strong in 
the metallic state (and actually modulated in the 3x3 
distortion), but counterproductive and thus weaker in 
the insulating state. Within the artificial -^73 geometry, 
results for Sn/Ge(lll) are on the whole similar to those 
in Sn/Si(lll), except that magnetism does not develop 
at U ~ 0, and the metal-insulator transition only occurs 
for U ~ A eV. All gaps and splittings are correspondingly 
smaller. 

We may crudely identify the magnetic 2D insulators 
near [/ ~ 4 eV with the actual MI states of the true 
surfaces, whose measured gaps are qualitatively similar. 



ranging from zero to hundreds of meV [17, |18|. Be- 
fore taking that identification seriously, we must how- 
ever compare the structural properties of the insulator 
with those of the competing phases, in particular with 
the 3x3 distorted metallic ones. To that end we re- 
peated all LSDA-f-U calculations in an enlarged 3x3 cell, 
where both distortive and ferromagnetic order parame- 
ters are allowed. The outcome for Sn/Si(lll) was un- 
eventful, and identically the same undistorted magnetic 
states were recovered for all U values. That confirms that 
a \/3 MI is indeed the LSDA-I-U predicted ground state 
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of Sn/Si(lll). Apart from magnetism, not yet investi- 
gated experimentally, this fully agrees with recent data 
by Modesti et al. [1^. 

The evolution is different in Sn/Ge(lll) where initially 
at ?7 ~ the ground state is metallic, nonmagnetic, but 
now 3x3 distorted, one Sn up by A/2, two down by 
A/2. The up-down amplitude A ~ 0.36 A and the energy 
gain ~ 9 meV/Sn are same as previously found long ago 



by similar methods [1, [13, U , 13] and hitherto believed 
to describe the true ground state of this surface. How- 
ever as shown in Figl2] (obtained by constrained struc- 
tural optimization along the distortion path) for increas- 
ing U the distortion A eventually disappears giving way 
to the same magnetic insulator previously found at ?7 ~ 4 
eV. Hence inclusion of correlations strongly modifies the 
ground state of Sn/Ge(lll) from a 3 x 3 distorted metal 
to a magnetic insulator with -^3 structural symmetry. 
This agrees with the STM and photoemission results of 
Cortes et al. below 20 K 
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In the MI phase, an unpaired electron is localized near 
each Sn adatom site. Ignoring anisotropy, these 1/2 spins 
form a 2D triangular Heisenberg-like magnet with in- 
tersite antiferromagnetic (AF) exchange coupling whose 
order of magnitude is set by J ~ t^/Ueff ~ 5 meV 
(for t - 0.05 eV, U^ff ~ 0.5 eV ^). Their ground 
state could thus be either a 120° noncollinear Neel state 
(irt or a spin liquid [13]. We made use of coUinear 
LSDA-I-U to calculate the actual J values for Sn/Si(lll) 
and Sn/Ge(lll). We carried out two separate 3x3 unit 
cell [/ = 4 eV calculations, one fully ferromagnetic (3/xb 
per cell), the other ferromagnetic (I^ib per cell), with 
two adatom spins up, one down [29|. The total energy 
change per 3x3 cell Ai? amounts to switching 6 bonds 
from antiferromagnetic to ferromagnetic, which for spins 
1/2 implies an energy loss of 9 J. We calculated Ai?=27 
and 42 meV, implying J ^ 3 and 5 meV for Sn/Si(lll) 
and Sn/Ge(lll) respectively. Because of the better ac- 
curacy of the ferro calculation, these are most likely un- 
derestimates, and we conclude that at least up to a hun- 
dred Kelvin the MI state of Sn/Si(lll) and Sn/Ge(lll) 
is likely to possess AF short range order, probably even 
outright Neel order, although renormalized by quantum 



fluctuations [27[. Hole spectral function calculations in 
the 2D triangular t- J model i30[] suggest that purely mag- 
netic 3x3 order will give rise to photoemission features 
typical of a 3 X 3 band folding, even in the absence of 
structural 3x3 periodicities. Fig[3] shows the LSDA+U 
electronic bands of Sn/Si(lll) ferrimagnetic state which 
exempli fy th at folding. The photoemission intensity near 
Ef at F [l3,[i9|] observed in Sn/Si(lll) is thus likely due 
to 3 X 3 AF order and consequent folding. The same 3x3 
magnetic folding should now be searched and detected 
at the true F-point of Sn/Ge(lll)'s low temperature MI 
state. 

These results are intriguing, and lead on to a number 
of interesting questions. 
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FIG. 3: (a) Sn/Si(lll): LDA-hU spin-up bands (solid line) 
and spin-down bands (dashed line) for the ferrimagnetic 3 x 
3 phase, (b) Total density of electronic states (solid line) 
showing the gap at the Fermi level (zero). 
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FIG. 4: Schematic phase diagram of Sn/Si(lll) in the 
temperature-interaction plane. The arrow describes the 
entropy-induced metallization due to AF order in the MI. 



First, why do the MI states evolve into metallic phases, 
beginning near 60 K and 20 K respectively in Sn/Si(lll) 
and in Sn/Ge(lll)? Moreover, why do these finite tem- 
perature metallic phases resemble the uncorrelated LDA 
ground states calculated at T=0? 

The driving force for the observed T-induced met- 
allization of MI states is, in analogy with V2O3 and 
2D fc - {ET)2Cu[N{CN)2]Cl [sij their lack of spin en- 
tropy, frozen out by AF short range order. If A_E ~ 
W{U /Uc — 1) represents a typical insulator-metal en- 
ergy difference per site (for a transition aX U — Uc), and 
AS{T) ~ the entropy-related free energy difference 
(assuming AF frozen spins and a charge gap in the MI), 
then the low-T Ml-metal phase boundary is predicted of 
the form: 
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which is sketched in FiglH and well brought out experi- 
mentally in 3D compounds. [3 Ij Here W ~ 0.3 eV is the 
bare bandwidth, and 7 is the electronic specific heat co- 
efficient (perhaps of order 0.01 J/mol K"^, or roughly 0.1 
meV/site K"^ as in 2D organics[3l|). We thus suggest 
that quasiparticle entropy drives these surfaces across 
the insulator - metal phase boundary, resurrecting the 
metallic phase, itself only metastable at T = 0. Four 
point conductance measurements [s^ could be of great 
help in ascertaining this scenario. 

The next and crucial if still speculative question is 
whether there is any chance to realize a 2D (power-law 
ordered) superconducting state in these surfaces. Q As 
indicated in FigHl in the general phase diagram of the 
triangular Hubbard lattice, as realized by organics under 
pressure[3lj - there is indeed a low-T d-wave supercon- 
ductor on the metallic side next to the MI phase. We thus 
propose that one should try to achieve superconductivity 
by metallization of Sn/Si(lll) below 60 K. Metallization 
could be attempted by e.g., reducing U, or by increasing 
the bandwidth, or by doping surface bands away from 
half-filling. The latter might be realized by alkali depo- 
sition; the former possibly through heavy doping of the 
Si bulk substrate. We note in passing that the d-wave 
superconducting state in a triangular lattice would prob- 
ably break rotational symmetry [33\ . an event readily 
observable by STM. 

One remaining unknown is the role of spin orbit cou- 
pling. While large in Sn, the largely Pz nature of the state 
will reduce its relevance. Possibly, some amount of mag- 
netocrystalline anisotropy will result. Depending on its 
sign, the spin 1/2 sites will turn from Heisenberg to either 
Ising (out of the surface plane; in this case, the ferrimag- 
netic state considered above actually corresponds to the 
true MI ground state) or XY (in the surface plane) . This 
aspect remains open for future investigation. 

In conclusion, in this Letter we provide a first theoret- 
ical background for the insulating ground state just ob- 
served in surface a-phases Sn/Si(lll) and Sn/Ge(lll), 
and obtain quantitative indications that the insulating 
phase observed in both should be of correlation ori- 
gin. LSDA+U calculations predict a stable magnetic and 
(Mott-Hubbard like) insulating \/3 phase in Sn/Si(lll), 
as well as the disappearance of the 3x3 distortion in 
Sn/Ge(lll), as is experimentally observed. Temperature 
induced metallization is argued to represent evidence for 
antiferromagnetism in the Mott state. Should metalliza- 
tion be provoked at sufficiently low temperatures, a 2D 
d-wave superconducting state could be achieved. A num- 
ber of experimental approaches, including surface doping, 
and four point conductance measurements are suggested 
for the future. 
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